The instability of black hole formation in gravitational collapse 
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We consider here the classic scenario given by Oppenheimer, Snyder, and Datt, for the gravitational collapse 
of a massive matter cloud, and examine its stability under the introduction of small tangential stresses. We show, 
by offering an explicit class of physically valid tangential stress perturbations, that an introduction of tangential 
pressure, however small, can qualitatively change the final fate of collapse from a black hole final state to a 
naked singularity. This shows instability of black hole formation in collapse and throws important light on the 
nature of cosmic censorship hypothesis and its possible formulations. The key effect of these perturbations is to 
alter the trapped surface formation pattern within the collapsing cloud and the apparent horizon structure. This 
allows the singularity to be visible, and implications are discussed. 
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The continual gravitational collapse of a massive matter 
cloud within the framework of general relativity was inves- 
tigated for the first time by the classic works of Oppenheimer 
and Snyder, and Datt (OSD) 1 1]. Such a treatment of dynami- 
cal collapse would be essential to determine the final fate of a 
massive collapsing star which shrinks catastrophically under 
the force of its own gravity when its internal nuclear fuel is 
exhausted. 

The outcome in the above case is seen to be a black hole de- 
veloping in the spacetime. As the gravitational collapse pro- 
gresses, an event horizon forms within the collapsing cloud 
and from the region within the horizon no material particles 
or light rays can escape, thus forming a black hole. The con- 
tinually collapsing star enters the horizon and finally ends up 
forming a spacetime singularity, which is hidden inside the 
black hole and which is unseen to all the outside observers in 
the universe. The matter and energy densities, spacetime cur- 
vatures, and all physical quantities blow up and take extreme 
values in the limit of approach to such a spacetime singularity. 

This classic picture became the foundation of an extensive 
theory and astrophysical applications of modern day black 
hole physics, further to the suggestion that all realistic mas- 
sive stars undergoing a continual gravitational collapse would 
have the same qualitative behaviour This means that, while 
the general theory of relativity necessarily implies the forma- 
tion of a spacetime singularity as the endstate for a massive 
collapsing star, such a singularity will always be necessarily 
hidden within a black hole. Such an assumption is known as 
the cosmic censorship hypothesis [2], and taking it to be valid, 
the theory and applications of black hole physics have devel- 
oped extensively in past many decades. 

The cosmic censorship has, however, remained an un- 
proved conjecture as yet in gravitation theory, despite numer- 
ous attempts to establish the same. Therefore, in past many 
years, much effort has also been devoted towards understand- 
ing and analyzing the final fate of a physically realistic dy- 
namical gravitational collapse scenario. The current status is. 
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despite much work in studying the censorship and its impli- 
cations, the issue of final fate of a complete gravitational col- 
lapse of a massive star remains far from being fully resolved. 
In particular, we need to formulate in a precise manner the 
conditions in gravitational collapse that would lead to the for- 
mation of black holes necessarily. We now know that under 
a wide variety of physically realistic situations, the collapse 
ends in a black hole or a naked singularity, depending on the 
initial conditions from which the collapse develops and the 
dynamical evolutions as allowed by the Einstein equations 
(see e.g. IH 01 and references therein). It is now clear that 
naked singularities are to be considered as a general feature of 
general relativistic physics and that they may develop as the 
end-state of collapse in a broad variety of physical collapse 
situations. 

It follows that a careful and extensive study of gravitational 
collapse phenomena in general relativity is the key to put the 
theory of black holes and their astrophysical implications on 
a firm footing. 

From such a perspective, we investigated here the effect of 
introducing small stress perturbations in the collapse dynam- 
ics of the classic Oppenheimer-Snyder-Datt gravitational col- 
lapse, an idealized model assuming zero pressure, which ter- 
minates in a black hole final fate. Our key purpose here is to 
study the stability of the OSD black hole under introduction 
of small tangential pressures. Clearly, stresses within a mas- 
sive collapsing star are very important physical forces to be 
taken into account while considering its dynamical evolution 
and the final fate of collapse (see for example |5]). 

We show here explicitly the existence of classes of stress 
perturbations such that the introduction of a smallest tangen- 
tial pressure within the collapsing OSD cloud changes the 
endstate of collapse to formation of a naked singularity, rather 
than a black hole. It follows that the OSD black hole is not 
stable under small stress perturbations within the collapsing 
cloud. As we point out below, this can also be viewed as per- 
turbing the spacetime metric of the cloud in a small way. Our 
work thus clarifies the role played by tangential stresses in a 
well known gravitational collapse scenario. The class of stress 
perturbations considered here, although specific, is physically 
reasonable and generic enough so as to provide a good insight 
into the stability of the OSD black hole. Clearly, such a result 
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provides an important insight into the structure of the censor- 
ship principle which as yet remains to be properly understood. 
This has also implications towards the physical consequences 
of final outcomes of a continual collapse, some of which are 
indicated in the concluding remarks. 

The general spherically symmetric line element describing 
the collapsing matter cloud can be written as. 



with the stress-energy tensor for a generic matter source given 
by, = -p, T; = pr, T« = = pe. The above is 
a general scenario, in that it involves no assumptions on the 
form of the matter or the equation of state. 

In order to decide on the stability or otherwise of the OSD 
model under the injection of small stress perturbations, we 
need to consider the dynamical development of the collaps- 
ing cloud, as governed by the Einstein equations. The visibil- 
ity or otherwise of the final singularity is determined by the 
behaviour of apparent horizon in the spacetime, which is the 
boundary of the trapped surface region that develops as the 
collapse progresses. First, we define a scaling function v{r, t) 
by the relation R — rv The Einstein equations for the 
above spacetime geometry can then be written as. 
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where the functions F and G are defined as, H — 

^-2v(r,v)-^2^ Q ^ g-2V(r,«)^/2^ jj^g ^^^^ 

tions in seven unknowns, namely p, p^, P6^ R, F, G, H. 
Here p is the mass-energy density, pr and po are the radial and 
tangential stresses respectively, R is the physical radius for the 
matter cloud, and F is the Misner-Sharp mass function. 

With the above definitions ofv,H and G, we can substitute 
the unknowns R, H with u, v. Without loss of generality, the 
scaling function v can be set v{ti,r) — 1 at the initial time 
ti = Q when the collapse commences. It then goes to zero 
at the spacetime singularity t^, which corresponds to i? = 0, 
i.e. we have v{ts, r) = 0. The above amounts to the scaling 
i? = r at the initial epoch, which is an allowed freedom. The 
collapse condition here is i? < throughout the evolution, 
which is equivalent to i; < 0. 

We can integrate (|4| by defining a suitably regular func- 
tion A{r,v) by = A^y{r,v)R' (the function A is defined 
in full generality here, while often the restriction to the class 
f = v{R), implying A{R) — v{R), is made, see e.g. |01). 
This gives G{r,t) = b{r)e^^^^^''"'. The arbitrary function of 
integration b{r) can be interpreted following the analogy with 
dust collapse models, where pressures vanish. It turns out to 
be related to the velocity of the collapsing shells, and once 
we write it as b{r) = 1 + r'^5o(r), we can see that values 



bo = const, in the dust limit correspond to the open (bo < 0), 
closed (bo > 0) or flat (bo = 0) Friedmann- Robertson- Walker 
models. The radial stress p,. and the energy density p are ob- 
tained from equations (|2]), once a specific choice for the mass 
function F{r, t) is made. The function v can be taken as the 
second free function for the system so that once a particular 
form of ly is specified, equation (O provides the tangential 
stress profile pe. Finally, from the equation of motion (|5]), 
we can integrate to obtain v{r, t), thus solving the system of 
Einstein equations. 

We can also invert the function v{r, t), which is monotoni- 
cally decreasing in t, to obtain the time needed by the matter 
shell at any radial value r to reach the event with a particular 
value V. We write the function t{r, v) from equation (|5]l as. 
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The time taken by the shell at r to reach the spacetime singu- 
larity at w = is then ts{r) = t{r, 0). 

Since t{r, v) is in general at least everywhere in the 
spacetime (because of the regularity of the functions in- 
volved), and is continuous at the center, we can write it as. 



t{r, v) = i(0, v) + rx{v) + 0{r^) 
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When t{r, v) is differentiable, we can make a Taylor expan- 
sion near the center r = 0. Here, t{0,v) is the above in- 
tegral evaluated at r = and xi"^) ~ ^L-o' 
point out below, the quantity x(0) plays an important role to- 
wards determining the nature of the final singularity of col- 
lapse. We consider collapse from a regular initial data, and 
so the Einstein equation (5) implies that the Misner-Sharp 
mass F(r, v) must go as near the center r = in order 
for the density to be regular at the center, and also to have 
t(0, v) well defined. Therefore, in general, F must have the 
form, F{r,v) = r'^M{r,v), where M is a suitably regular 
function. Then, by continuity, the time for the shell located 
at any r close to center to reach the singularity is given as, 
tg{r) = ts{0) + rx(0) + 0{r^). Basically, this means that 
the singularity curve should have a well-defined tangent at the 
center. Regularity at the center also implies that the metric 
function v cannot have constant or linear terms in r in a close 



neighborhood of r = 0, and it must go as v 



near the 



center. Therefore the most general choice of the free function 



1^ IS, 



i^(r, v) = r'^g{r, v) 



(8) 



Since g{r, v) is a regular function (at least C^), it can be writ- 
ten near r = as. 



g{r, v) = goiv) + gi{v)r + g2{v)r'^ 



(9) 



We would now like to investigate how the OSD gravita- 
tional collapse scenario, which is a homogeneous pressureless 
dust cloud collapsing to give rise to a black hole, gets altered 
when small stress perturbations are introduced in the dynami- 
cal evolution of collapse. 
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To that end, we first note that the dust scenario is obtained 
if Pr — pe = in the above. In that case, from equation 
(|3]l it follows that = and that together with the condition 
j^(0) = gives 1^ = identically. These models have been 
widely studied in the literature, and it is seen that for generic 
dust collapse the final outcome can be either a black hole or a 
naked singularity, depending on the nature of the initial den- 
sity and velocity profiles of the collapsing matter shells |8]. In 
the OSD collapse to a black hole, the trapped surfaces or the 
apparent horizon in the spacetime develop much earlier before 
the formation of the final singularity of collapse. On the other 
hand, when inhomogeneities are allowed in the initial density 
profile, such as a higher density at the center of the star, then 
the trapped surface formation is delayed in a natural manner 
within the collapsing cloud and the final singularity becomes 
visible to faraway observers in the universe [gj. 

The OSD case is obtained from above when we further as- 
sume that the collapsing dust is necessarily homogeneous at 
all epochs of collapse. This is of course an idealized scenario 
because realistic stars would have typically higher densities at 
the center, which slowly falls off with increasing radius, and 
they also would have non-zero internal stresses. Specifically, 
the conditions that must be imposed to obtain the OSD case 
from the above are given by, 

(a) M = Mo , 

(b) V ^ v(t) , 

(c) boir) = k . 

Then we have F' = 3Mor^, R' ~ v, and the energy 
density is homogeneous throughout the collapse with p — 
p{t) — iMo/v^. The spacetime geometry then becomes the 
Oppenheimer-Snyder metric. 



1 + kr^ 
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where the function v{t) is solution of the equation of motion, 
^ = -^Z (Mo/v) + k, obtained from Einstein equation {5^. 
In this case we get x(0) = identically. All the matter shells 
then collapse into a simultaneous singularity (due to condition 
(b)), which is necessarily covered by the event horizon that 
developed in the spacetime at an earlier time, thus giving rise 
to a black hole. 

To examine the effect of introducing stress perturbations 
in the above scenario and to study the models thus obtained 
which are close to the Oppenheimer-Snyder in this sense, we 
need to relax and perturb one or more of the above conditions 
(a), (b) or (c). 

If the collapse outcome would not to be a black hole, the 
final singularity of collapse cannot be simultaneous. We are 
thus led to relax condition (b) above, allowing v — v{t, r), 
rather than v = v{t) only. At the same time, in order not to 
depart too much from the OSD model, we keep (a) and (c) 
unchanged. This also brings out more clearly the role played 
by the stress perturbations in the model. 

In terms of the spacetime metric ([T]l, while the metric func- 
tion ^{t, r) must be identically vanishing for the dust case, the 



above amounts to allowing for small perturbations in i^, and 
allowing it to be non-zero now. This is equivalent to intro- 
ducing small stress perturbations in the model, and we show 
below how that affects the apparent horizon developing in the 
collapsing cloud. 

We note immediately that taking M — Mq leads to = 
r^Mo. We have R' = v + rv' — !• v for r — ^ and therefore 
we get A „ = v' /v. With the expansion near r — Q for both 
A and g we get the relation between the coefficients of the 
expansion of g and those for the expansion of A. Integrating 
© in the small r limit we thus obtain G(r, t) = fe(r)e^''('''"^. 
The radial stress pr vanishes in this case as F = 0, while the 
tangential pressure, obtained from equation (O, has the form, 
Pe — Pir^ + P2i'^ + where pi,p2 are naturally evaluated 
in terms of coefficients of to, g, and R and its derivatives. 
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Here the choice of sign of the functions go and gi is enough 
to ensure positivity or negativity of pg. 

We note that scenarios with vanishing radial stresses but 
non-vanishing tangential stresses have been considered in 
past, with the most physically significant model (though not 
the o nly relevant one) being the so called 'Einstein cluster' 
(see lUOtl '). which describes a cloud of collapsing counter ro- 
tating particles. Naked singularities and black holes are found 
to arise as the endstate of such models, depending on the ini- 
tial density, velocity and stress configurations 1 11]. 

The first order coefficient x in equation of the time curve of 
the singularity tg (r) is now obtained as 



X(0) = - 
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As mentioned above, it is x(0) that governs the nature of the 
singularity curve, and whether it is increasing or decreasing 
away from the center Clearly, it is the matter initial data in 
terms of density and stress profiles, the velocity of the collaps- 
ing shells, and the allowed dynamical evolutions that govern 
and fix the value of x(0). 

The quantity x(0) also governs the behaviour of apparent 
horizon and trapped surface formation, as we show below, 
which in turn governs the nakedness or otherwise of the sin- 
gularity. The equation for the apparent horizon is given by 
F / R = 1. It is analogous to that of the dust case since 
F/R = rM/v in both cases |9]. So the apparent horizon 
curve rah{t) is given by rj,^ = f^, with Vah = v{rah{t),t), 
which can also be inverted as a time curve tahif)- The visi- 
bility of the singularity at the center of the collapsing cloud to 
faraway observers is determined by the nature of this apparent 
horizon curve which is given by. 



tah{r) = ts{r) - 



zdv 



Mo 



(13) 



where ts(r) is the singularity time curve, whose initial point 
is to = ts{0)- Near r = the equation (fT3b becomes. 



tahir) ^ to + xiO)r + oir^) 



(14) 
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From the above, it is now easy to see how the stress pertur- 
bation affects the time of formation of the apparent horizon, 
and therefore the formation of a black hole or naked singular- 
ity. A naked singularity typically occurs as collapse endstate 
when a co-moving observer at fixed r does not encounter any 
trapped surfaces till the time of singularity formation. For a 
black hole to form, trapped surfaces develop before the singu- 
larity, so it is needed that, 

tahir) < to for r > 0, near r = . (15) 

It is clear that for all functions gi{v) for which x(0) is pos- 
itive, this condition is violated and the apparent horizon is 
forced to appear after the formation of the central singular- 
ity. The apparent horizon curve then initiates at the central 
singularity r — at t — to and increases with increasing r, 
moving to the future, i.e. tah > to for r > near the cen- 
ter The behaviour of outgoing families of null geodesies has 
been analyzed in detail in such a case when x(0) > and we 
know that geodesies terminate at the singularity in the past. 
Thus timelike and null geodesies come out from the singular- 
ity, making it visible to external observers llT2ll . 

It follows that gi is the term in the stresses pg which de- 
cides the black hole or naked singularity final fate. We can 
choose it to be arbitrarily small, and we now see how intro- 
ducing a generic tangential stress perturbation in the model 
would change drastically the final outcome of collapse. For 
all non-vanishing tangential stresses with go = and gi < 0, 
even the slightest perturbation of the Oppenheimer-Snyder- 
Datt scenario, injecting a small tangential stress would result 
in a naked singularity. The space of all functions gi that make 
x(0) positive, which includes all the strictly negative func- 
tions gi, causes the collapse to end in a naked singularity. We 
note that while this is an explicit example, by no means this is 
the only class. 

The remarkable feature of this class is that it corresponds to 
a collapse model for a simple and straightforward perturbation 
of the Oppenheimer-Snyder-Datt spacetime metric, where the 
geometry near the center can be written as, 

ds^ = -(1 - 2g,r')dt^ + _(^±I!Ql_dr2 + r^'dQ^ , 

1 + Kr^ — 2gir'^ 

(16) 

The metric above satisfies Einstein equations in the neigh- 
borhood of the center of the cloud when the function gi{v) 
is small and bounded. We could take for example, < 
|(7i(ii)| < e, so that the smaller we take the parameter e the 
bigger will be the radius where the approximation is valid. 
The function v{r,t) above is governed by the equation of 
motion (|5]i which in the small r limit becomes, dv/dt = 
(1 - gi{v)r^){^ + k- 25i(i;)r)i/2. Finally, x(0) in this 
case is given by equation (fTZt with go = 0, and in certain 
cases can also be integrated. 

We note that any realistic matter model must satisfy some 
energy conditions ensuring the positivity of mass and energy 
density. In general, the weak energy condition implies restric- 
tions on the density and pressure profiles. The energy density 



as given by the second of equations (O must be positive. Since 
R is positive, to ensure positivity of p we require F > and 
R' > 0. The choice of positive M(r) (which obviously holds 
for Mo > and is physically reasonable) ensures positivity 
of the mass function. Then i?' > is a sufficient condition 
for the avoidance of shell crossing singularities. The tangen- 
tial stress can be written from ((Sj where pr = 0, and is given 
by Pe — \^P^' ■ So the sign of the function v' determines 
the sign of pg. Positivity of p + pg is then ensured for small 
values of r throughout collapse for any form of pg. In fact, 
regardless of the values taken by M and g, there will always 
be a neighborhood of r = for which \pg\ < p and therefore 
P + Pe > 0. 

The black hole and naked singularity outcomes of gravita- 
tional collapse are very different from each other physically, 
and would have quite different observational signatures. In the 
naked singularity case we have the possibility to observe the 
physical effects happening in the vicinity of the ultra dense 
regions that form in the very final stages of collapse. How- 
ever, in a black hole scenario, such regions are necessarily hid- 
den within the event horizon. The fact that a slightest stress 
perturbation of the OSD collapse could change the outcome 
drastically, taking it from a black hole to naked singularity 
formation, means that the naked singularity final state for a 
collapsing star must be studied carefully to deduce its physi- 
cal consequences which are not well understood so far. 

The existence of subspaces of collapse solutions as we have 
shown here, that go to a naked singularity final state rather 
than a black hole, in the arbitrary vicinity of the OSD black 
hole, presents an intriguing scenario. It gives an idea of the 
richness of the structure present in gravitation theory and the 
complex solution space of Einstein equations which are a 
complicated set of highly non-linear partial differential equa- 
tions. What we see here is there are classes of stress pertur- 
bations such that an arbitrarily small change from the OSD 
model is a solution going to naked singularity. In this sense, 
this manifests an instability in the black hole formation pro- 
cess in gravitational collapse. This also provides an intriguing 
insight into the nature of cosmic censorship, namely that the 
collapse must be properly fine-tuned necessarily if it is to pro- 
duce a black hole only as the final endstate. 

Traditionally it was believed that the presence of stresses 
or pressures in the collapsing matter cloud would increase the 
chance of black hole formation, thereby ruling out dust mod- 
els that were found to lead to a naked singularity as collapse 
endstate. That is no longer the case. The model described 
here not only provides a new class of collapses ending in a 
naked singularity, but more importantly, shows how the bi- 
furcation line that separates the phase space of 'black hole 
formation' from that of the 'naked singularity formation' runs 
directly over the simplest and most studied of black hole sce- 
narios such as the OSD model, thus making it unstable under 
perturbations. 
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